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INTRODUCTION 

Metastasis  is  one  of  the  biggest  problems  in  cancer  therapy.  E-cadherin  (CDH1)  is 
known  to  be  involved  in  epithelial-mesenchymal  transition  (EMT)  process,  thus  playing 
an  important  role  during  the  tumor  metastatic  progression  [1-3].  It  is  reported  that  losing 
CDH1  expression  in  conjunction  with  EMT  might  contribute  to  the  development  of 
multidrug  resistance  (MDR)  to  chemotherapy  [4-6],  Thus,  understanding  the  CDH1 
down-regulation  mediated  EMT  and  its  relationship  with  MDR  will  have  great  benefit  for 
improving  chemotherapy  efficiency.  Most  of  the  present  in  vitro  research  analyzing  the 
functions  of  CDH1  and  its  involvement  in  mediating  metastasis  have  been  carried  out  in 
2D  cultured  cell  lines.  However,  it  has  been  reported  that  cells  in  2D-culture  systems 
have  different  gene  expression  patterns  than  in  vivo  cells  [7-11].  By  comparison,  the  3D 
culture  systems  provide  more  accurate  molecular  pictures  of  human  disease.  The 
establishment  of  3D-culture  systems  will  be  a  better  way  to  mimic  the  tissue  samples  and 
provide  us  with  more  information  that  recapitulates  human  tissues. 

In  this  study,  we  aim  to  establish  a  3D  model  that  can  simulate  the  in  vivo  cancer 
models  and  can  be  induced  to  down  regulate  its  CDH1  expression  when  stimulated,  and 
thus  can  be  used  for  simulation  of  EMT  processes.  We  selected  CDH1  as  a  target  for 
manipulation  because  its  expression  is  known  to  cause  pronounced  chemical  and 
phenotypic  changes  in  epithelial  cells.  The  CDH1  knock  down  needs  to  be  inducible  as 
expression  of  CDH1  is  essential  during  formation  of  the  3D  structures  and  the  CDH1  will 
be  knocked  down  after  the  3D  structures  are  fully  formed.  We  hypothesize  that  by 
reducing  the  expression  of  CDH1  in  colorectal  cancer  cell  lines  cultured  in  3D-structures, 
the  colorectal  cancer  cells  will  display  EMT  characteristics  and  gain  drug-resistance 
compared  to  CDH1  (+)  colorectal  cancer  cells.  We  hypothesize  that  CDH1  reduction  will 
induce  changes  in  downstream  pathways,  including  proteins  in  pathways  of  cell  survival 
(PI3K,  GSK3b,  EGFR)  apoptosis  (caspases  family  proteins,  Bcl-2  family  proteins),  and 
cell  migration  (Rho  family  proteins  like  Cdc42,  Racl,  RhoA,  WASP  family  proteins, 
and  PAR/aPKC  complexes).  Our  research  will  be  the  first  to  connect  EMT  to  drug 
resistance  in  3D-culture  systems. 

During  the  first  year  of  this  research,  we  put  our  effort  mainly  into  establishing 
and  optimizing  the  inducible  CDH1  knock  down  model  with  3D  culturing  system,  and 
evaluating  its  characteristics  and  performance.  During  the  second  year  of  this  research, 
we  put  our  effort  mainly  into  establishing  and  optimizing  the  methodology  for  proteomic 
and  phosphoproteomic  analysis,  and  started  to  collect  data  for  evaluating  the  functions  of 
CDH1  for  drug  resistance  with  mass  spectrometer.  These  data  will  be  analyzed  for 
changes  in  the  signal  transduction  pathways  mentioned  above. 
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ACCOMPLISHMENTS 

1.  (Task  2a)  Establish  the  SILAC  method.  (7/1/2013-8/30/2013) 

In  our  proposed  research,  we  will  perform  quantitative  experiments  to  study  the 
molecular  changes  of  3D  structures  in  response  to  CDH1-KD  induced  drug  resistance. 
We  will  use  the  Stable  Isotope  Labeling  of  Amino  Acids  in  Cell  Culture  (SILAC) 
strategy  to  label  3D  structures  under  different  conditions,  and  then  mix  the  samples 
together  and  perform  LC-MS/MS  analysis.  To  prepare  for  the  quantitative  analysis,  we 
first  analyzed  the  labeling  efficiency  of  the  SILAC  method.  SILAC  labels  were  checked 
both  in  2D  and  3D  culturing  systems.  The  results  show  that  in  both  2D  and  3D  culturing 
systems,  for  the  cells  cultured  in  media  supplemented  with  heavy  amino  acids,  around 
98%  of  peptides  were  labeled  with  heavy  amino  acids.  Lor  the  cells  cultured  in  media 
supplemented  with  light  amino  acids,  more  than  99.8%  of  the  peptides  have  light  amino 
acids,  and  the  presence  of  heavy  amino  acids  is  less  than  0.2%  (Figure  1A).  The  protein 
amounts  in  each  sample  were  determined  by  BCA  assay.  However,  because  the  BCA 
assay  can  sometimes  generate  some  error,  to  ensure  that  in  the  real  experiments  the  light 
and  heavy  samples  can  be  mixed  with  exactly  a  1:1  ratio,  a  pre-experiment  was 
performed  to  determine  the  actual  heavy-to-light  mixing  ratio  (vol/vol).  The  heavy- 
labeled  3D  structures  and  light-labeled  2D  layers  were  mixed  with  different  volume  ratios 
to  make  a  1 : 1  total  protein  ratio  determined  by  BCA  assay  (Figure  IB).  Six  samples  were 
generated  with  different  vol/vol  ratios  and  tested  with  LC-MS/MS  with  60  min  run  (blue 
line  in  Figure  IB)  and  120  min  run  (red  line  in  Figure  IB).  The  heavy-to-light  ratio 
determined  by  LC-MS/MS  showed  a  linear  correlation  to  the  vol/vol  ratios  determined  by 
BCA  assay.  The  light-labeled  3D  structures  and  heavy-labeled  2D  structures  were  tested 
with  the  same  strategy  and  showed  similar  results  (Figure  1C).  From  these  linear 
correlations  we  determined  the  actual  vol/vol  ratio  to  verify  that  the  peptide  heavy-to- 
light  ratio  was  1:1  for  both  60  min  and  120  min  runs.  We  will  use  the  average  of  these 
two  numbers  to  make  the  heavy  and  light  sample  mixtures. 
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Figure  1 
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Figure  1.  SILAC  labeling  efficiency.  (A)  SILAC  heavy  and  light  amino  acid  label 
efficiencies  are  shown  in  bar  graph.  (B,  C)  Validation  of  the  ratio  (vol/vol)  to  mix 
heavy  and  light  samples  to  make  1:1  SILAC  mixtures.  (B)  is  for  the  mixture  of 
heavy-labeled  3D  structures  and  light-labeled  2D  layers,  and  (C)  is  for  the  mixture 
of  light-labeled  3D  structures  and  heavy-labeled  2D  layers. 
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2.  (Task  2a  -  Extension)  Optimizing  phosphoproteomic  methods  for  phosphoprotein 
quantification.  (9/1/2013-1/31/2014) 

Phosphorylation  and  dephosphorylation  of  key  cellular  factors  and  the 
deregulation  of  kinase  signaling  pathways  are  commonly  associated  with  various  cancers. 
Among  the  many  biochemical  mechanisms  involved  in  cellular  signaling,  protein 
phosphorylation  is  one  of  the  most  common  control  mechanisms  to  regulate  a  variety  of 
biological  processes,  including  cell  proliferation,  migration,  DNA  reparation,  and 
apoptosis.  Over  the  past  decades,  great  efforts  were  made  to  estimate  the  phosphorylation 
level  of  signaling  pathways  in  different  types  of  cancers.  Phosphorylation  events  were 
reported  in  lung  cancer,  breast  cancer,  colorectal  cancer,  pancreatic  cancer,  and  prostate 
cancer.  We  published  a  manuscript  in  Journal  of  Proteome  Research  in  the  first  year  of 
the  project  introducing  a  new  phosphoproteomic  enrichment  workflow  that  works  better 
than  the  gold  standard.  We  are  further  improving  that  method  with  the  following  study. 

In  order  to  make  it  possible  to  identify  as  many  phosphoproteins  and  phosphosites 
as  possible,  we  also  optimized  the  methods  for  mass  spectrometric-based 
phosphoproteomic  studies.  In  our  experiments,  we  compared  two  materials  including 
IMAC  and  TiCh  that  are  widely  used  to  enrich  phosphopeptides  from  the  peptide 
mixtures  for  phosphoproteomic  studies. 

We  first  optimized  the  TiCh  or  IMAC  bead  to  peptide  ratios  to  achieve  the  best 
performance  of  the  beads  (Figure  2).  The  results  show  that  TiC>2  beads  perform  best  with 
the  TiC>2  to  peptide  ratio  ranging  from  2:1  to  8:1,  which  can  achieve  more  than  80% 
phosphopeptide  enrichment  efficiency  without  losing  total  phosphopeptide  identifications 
(Figure  2A).  The  IMAC  beads  show  the  best  performance  when  the  peptide  to  IMAC 
ratio  is  100:1,  which  shows  the  highest  enrichment  efficiency  of  more  than  60%  (Figure 
2B). 
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Figure  2 
A 


Optimization  ofTi02  to  Peptide  Ratio 
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Peptide:IMAC  (mg/ml) 

i  i Total  Peptides  i  i  Phosphopeptides  — ♦-*  Phosphopeptide  Ratio  (%) 

Fig  2.  Optimization  of  the  peptide  to  bead  ratio  for  phosphopeptide  enrichment 
from  whole  cell  lysate.  (A)  Optimization  of  Ti02  to  peptide  ratio.  (B). 
Optimization  of  peptide  to  IMAC  ratio.  Identified  total  peptide  numbers  with 
different  Ti02-to-peptide  ratios  (u g : u g )  or  peptide-to-IMAC  ratios  (ug:ul)  are 
shown  in  the  bar  graph.  Percentages  of  identified  phosphopeptides  with  each 
condition  are  shown  with  the  dashed  line. 


With  the  optimized  material  to  peptide  ratios,  we  performed  direct  enrichment  of 
phosphopeptides  with  IMAC  or  Ti02  beads  from  whole  cell  lysates.  As  starting  material, 
3  mg  of  MCF-10A  cell  lysate  were  used  for  phosphopeptide  enrichment.  According  to 
the  optimized  material  to  peptide  ratios,  we  added  30  pi  of  the  IMAC  or  6  mg  of  TiCf 
bead  slurry  to  the  peptide  mixture  to  perform  phosphopeptide  enrichment.  The  results 
showed  that  a  total  of  4811  peptides  were  identified  with  the  first  round  of  Ti02 
enrichment  (TiCh-l)  and  among  which  3918  were  phosphopeptides,  with  a 
phosphopeptide  enrichment  ratio  of  81.4%  (Figure  3B&C).  A  total  of  4661  peptides 
were  identified  with  the  first  round  of  IMAC  enrichment  (IMAC-1),  and  among  these 
peptides,  3699  were  phosphopeptides,  with  a  phosphopeptide  enrichment  ratio  of  79.4% 
(Figure  3D&E).  PhosphoRS  was  used  to  evaluate  the  accuracy  of  phosphosites 
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localization  [12].  As  reported  previously,  phosphopeptides  with  PhosphoRS  score  >  99% 
are  considered  to  have  confident  phosphosite  localization  [12].  We  only  counted 
phosphopeptides  with  confident  localization  for  mono-  or  multi-phosphopeptide 
identification,  as  well  as  for  calculating  total  number  of  phosphosites  identified  in  all 
species.  TiC^-l  enrichment  generated  2631  phosphopeptides  with  confident  localization 
information,  among  which  2382  were  mono-phosphopeptides,  while  249  were  multi- 
phosphopeptides  (Figure  3B&C).  IMAC-1  enrichment  generated  2608  phosphopeptides 
with  confident  localization  information,  among  which  2166  were  mono-phosphopeptides, 
and  442  were  multi-phosphopeptides  (Figure  3D&E). 


Figure  3 


Flowthrough  Flowthrough 


Fig  3.  Distribution  of  phosphopeptides  identified  by  multi-step  Ti02  or  IMAC 
enrichment.  (A)  Experiment  workflow  for  multiple  rounds  of  enrichment.  (B) 
Numbers  of  identified  phosphopeptides  in  the  first,  second,  and  third  round  of  Ti02 
enrichment  (phosphoRS  score  >  99%)  with  1,  2,  3,  or  4  phosphate  groups  are 
shown  in  the  bar  graph  with  light  purple,  blue,  yellow  and  red  colors.  Unique 
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phosphosite  numbers  are  calculated  from  phosphopeptides  with  phosphoRS  score  > 
99%  and  are  shown  with  the  dashed  line.  (C)  Numbers  of  identified 
phosphopeptides  in  the  first,  second,  and  third  round  of  Ti02  enrichment  with 
phospho-serine  (blue),  phospho-threonine  (green),  or  phospho-tyrosine  (red)  sites, 
and  those  without  confident  localization  (phosphoRS  score  <  99%,  dark  purple)  are 
shown  in  the  bar  graph.  The  percentages  of  identified  phosphopeptides  from  the 
total  identified  peptides  in  each  round  of  IMAC  enrichment  are  shown  with  the 
dashed  line.  (D)  Numbers  of  identified  phosphopeptides  in  the  first,  second,  and 
third  round  of  IMAC  enrichment  (phosphoRS  score  >  99%)  with  1,  2,  3,  or  4 
phosphate  groups  are  shown  in  the  bar  graph  with  light  purple,  blue,  yellow  and  red 
colors.  Unique  phosphosite  numbers  are  calculated  from  phosphopeptides  with 
phosphoRS  score  >  99%  and  are  shown  with  the  dashed  line.  (E)  Numbers  of 
identified  phosphopeptides  in  the  first,  second,  and  third  round  of  IMAC 
enrichment  with  phospho-serine  (blue),  phospho-threonine  (green),  or  phospho- 
tyrosine  (red)  sites,  and  those  without  confident  localization  (phosphoRS  score  < 
99%,  dark  purple)  are  shown  in  the  bar  graph.  The  percentages  of  identified 
phosphopeptides  from  the  total  identified  peptides  in  each  round  of  IMAC 
enrichment  are  shown  with  the  dashed  line. 


To  evaluate  how  complete  a  single  enrichment  would  be  at  extracting 
phosphopeptides  from  the  whole  cell  lysate,  we  then  carried  out  a  second  round  of 
enrichment  from  the  flow-through  of  the  first  round  of  enrichment.  We  repeated  the 
process  for  a  third  round,  with  enrichment  from  the  flow-through  of  the  second  round  of 
enrichment.  In  the  second  round  of  TiC>2  enrichment,  1899  total  phosphopeptides  were 
identified.  Of  these,  1266  phosphopeptides  were  also  identified  with  the  first  round  of 
Ti02  enrichment,  while  525  phosphopeptides  were  unique  to  the  second  round 
enrichment.  With  the  third  round  of  Ti02  enrichment,  a  total  number  of  538 
phosphopeptides  were  identified,  among  which  350  phosphopeptides  had  been  identified 
with  the  first  round  of  Ti02  enrichment  and  430  phosphopeptides  overlapped  with  the 
second  round  of  TiC>2  enrichment.  Only  an  additional  80  unique  phosphopeptides  were 
detected  with  the  third  round  of  Ti02  enrichment  (Figure  4A).  In  the  second  round  of 
IMAC  enrichment,  3079  total  phosphopeptides  were  identified.  Of  these,  1791 
phosphopeptides  were  also  identified  with  the  first  round  of  IMAC  enrichment,  while 
1288  phosphopeptides  were  unique  to  the  second  round  phosphopeptides.  With  the  third 
round  of  IMAC  enrichment,  a  total  number  of  1531  phosphopeptides  were  identified, 
among  which  947  phosphopeptides  had  been  identified  with  the  first  round  of  IMAC 
enrichment  and  1252  phosphopeptides  overlapped  with  the  second  round  of  IMAC 
enrichment.  Only  an  additional  223  unique  phosphopeptides  were  detected  with  the  third 
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round  of  IMAC  enrichment  (Figure  4B).  All  the  above  results  indicate  that  three  cycles 
of  Ti02  or  IMAC  enrichment  would  be  enough  to  enrich  most  of  the  phosphopeptides 
from  the  whole  cell  lysates. 


Figure  4 


Fig  4.  Venn  diagram  of  phosphopeptides  identified  by  multi-step  Ti02  or  IMAC 
enrichment.  (A)  Overlaps  of  the  nonredundant  phosphopeptides  identified  in  each 
round  of  Ti02  enrichment  are  shown  in  the  Venn  diagram.  (D)  Overlaps  of  the 
nonredundant  phosphopeptides  identified  in  each  round  of  IMAC  enrichment  are 
shown  in  the  Venn  diagram. 

We  next  compared  the  differences  of  phosphopeptides  that  are  enriched 
with  the  multiple  rounds  of  Ti02  or  IMAC.  A  total  number  of  4727  unique 
phosphopeptides  were  identified  with  3  rounds  of  IMAC  enrichment,  while  a  total 
number  of  4537  unique  phosphopeptides  were  identified  with  3  rounds  of  TiCh 
enrichment,  showing  that  both  IMAC  and  TiC>2  have  the  equivalent  ability  to  enrich 
phosphopeptides  from  whole  cell  lysate  (Figure  5A).  IMAC  enriched  2793  mono- 
phosphopeptides  and  428  multi-phosphopeptides,  while  TiCh  enriched  2722  mono- 
phosphopeptides  and  259  multi-phosphopeptides,  indicating  that  IMAC  is  more  efficient 
for  enriching  multi-phosphopeptides.  The  distribution  of  phosphor-Ser,  phosphor-Thr  and 
phosphor-Tyr  with  IMAC  or  TiC>2  enrichment  is  shown  in  Figure  5C,  and  does  not  show 
much  difference  with  these  two  different  enrichment  methods.  However,  the  overlap 
between  phosphopeptides  identified  with  IMAC  or  TiC>2  shows  only  about  1/3  of  the 
phosphopeptides  were  identified  with  both  IMAC  and  TiC>2  enrichment  (Figure  5D). 
Among  the  overall  identified  phosphopeptides,  2364  were  unique  to  IMAC  enrichment, 
2174  were  unique  to  TiC>2  enrichment,  and  2364  were  identified  with  both  methods 
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(Figure  5D).  In  order  to  further  investigate  how  the  IMAC  enrichment  is  different  from 
Ti02  enrichment,  we  divided  the  overall  identified  phosphopeptides  into  3  subsets,  the 
first  subset  consists  of  phosphopeptides  that  are  uniquely  identified  by  IMAC  enrichment, 
the  second  subset  contains  phosphopeptides  that  are  uniquely  identified  by  Ti02 
enrichment,  and  the  third  subset  is  phosphopeptides  that  are  identified  both  with  IMAC 
and  Ti02  enrichment.  The  results  clearly  show  that  the  phosphopeptides  that  are  uniquely 
identified  with  IMAC  enrichment  contain  3  times  as  many  as  multi-phosphopeptide 
compared  with  those  that  are  uniquely  identified  with  HO2  (Figure  5E).  In  comparison, 
the  percentage  of  phosphor-Ser,  phosphor-Thr,  and  phosphor-Tyr  does  not  show  much 
difference  among  these  3  parts  (Figure  5F).  This  result  is  significant  as  multi- 
phosphorylated  peptides  are  biologically  important  but  obtain  hard  to  detect  and 
characterize.  With  this  approach,  we  have  determined  an  optimized  approach  to  identify 
multi-phosphorylated  peptides  from  complex  biological  samples. 


Figure  5 
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Figure  5 
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Fig  5.  Comparison  between  Ti02  and  IMAC  enriched  phosphopeptides.  (A) 
Distribution  of  non-phosphopeptides  and  phosphopeptides  identified  with  Ti02  or 
IMAC  enrichment.  (B)  Distribution  of  phosphopeptides  with  1,  2,  3,  or  4  phosphate 
groups  identified  with  Ti02  or  IMAC  enrichment.  (C)  Distribution  of 
phosphopeptides  with  phospho-serine,  phospho-threonine,  or  phospho-tyrosine 
sites  identified  with  Ti02  or  IMAC  enrichment.  (D)  Overlaps  of  the  nonredundant 
phosphopeptides  identified  with  Ti02  or  IMAC  enrichment  are  shown  in  the  Venn 
diagram.  (E)  Distribution  of  single  or  multi-phosphopeptides  that  are  identified  as 
IMAC  unique,  Ti02  unique,  or  with  both  enrichment  methods.  (F)  Distribution  of 
phosphopeptides  with  phospho-serine,  phospho-threonine,  or  phospho-tyrosine 
sites  that  are  identified  as  IMAC  unique,  Ti02  unique,  or  with  both  enrichment 
methods. 


We  next  evaluated  whether  the  differences  of  phosphopeptides  enriched  by  IMAC 
or  Ti02  result  from  the  enrichment  step,  or  from  the  elusion  step.  We  used  2  rounds  of 
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Ti02  to  further  enrich  phosphopeptides  from  the  flow-through  of  3  rounds  of  IMAC 
enrichment,  and  we  used  2  rounds  of  IMAC  to  further  enrich  phosphopeptides  from  the 
flow-through  of  3  rounds  of  Ti02  enrichment  (Figure  6A).  The  results  show  that  with 
further  TiC>2  enrichment  from  the  flow-through  of  IMAC  enrichment,  the  identified 
phosphopeptides  doesn’t  show  much  overlap  with  those  that  are  identified  with  TiC>2 
enrichment  from  the  whole  cell  lysate  (Figure  6B-E),  indicating  that  the  loss  of 
phosphopeptide  identification  is  not  due  to  the  incomplete  enrichment  from  the  cell  lysate, 
but  is  due  to  the  incomplete  elution  from  IMAC  or  HO2  materials. 


Figure  6 
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Figure  6 


C  T123  Unique  1311 12 


2000 

1800 

1600 

1400 

1200 

1000 

800 

600 

400 

200 


0 

I123T12 

Total 

Overlap 
with  1123 
Unique 

Overlap 
with  T123 
Unique 

Overlap 
with  IT 
Overlap 

Unique  in 
I123T1T2 

|t123I  112 

1888 

295 

265 

1002 

326 

Fig  6.  Distribution  of  phosphopeptides  identified  by  multi-step  Ti02  following 
multi-step  IMAC  enrichment.  (A)  Experiment  workflow  for  multiple  rounds  of 
enrichment.  (B-D)  Overlap  of  Ti02  enriched  phosphopeptides  from  IMAC  flow¬ 
through  with  phosphopeptides  in  IMAC  unique,  Ti02  unique,  or  IMAC-Ti02 
overlap  parts.  (E)  Distribution  of  phosphopeptides  identified  in  each  round  of  Ti02 
or  IMAC  enrichment. 


To  evaluate  the  contribution  of  each  round  of  enrichment  to  the  total  unique 
phosphopeptide  identification,  for  each  round  of  enrichment,  the  total  number  of 
identified  phosphopeptides  and  the  number  of  unique  phosphopeptides  are  summarized  in 
Figure  7A.  The  results  show  that  the  first  round  of  IMAC  or  Ti02  enrichment  gave  the 
highest  number  of  both  total  phosphopeptide  and  unique  phosphopeptide  identification, 
while  the  identification  decreases  for  the  second  and  third  round  of  enrichment  with  the 
same  material  (Figure  7A).  However,  by  switching  the  enrichment  material  to  the  other 
media,  the  identification  numbers  increased  again  (Figure  7A).  By  integrating  IMAC 
with  Ti02  enrichment,  the  total  number  of  identified  phosphopeptides  has  reached  8747 
with  just  2  days  of  workload.  The  contribution  of  each  round  of  enrichment  for 
identification  of  unique  phosphopeptides  is  summarized  in  Figure  7B. 
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Figure  7 
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Fig  7.  Distribution  of  phosphopeptides  identified  in  each  round  of  Ti02  or  IMAC 
enrichment.  (A)  Distribution  of  total  phosphopeptides  (blue)  and  unique 
phosphopeptides  (pink)  identified  in  each  round  of  Ti02  or  IMAC  enrichment.  (B) 
Contribution  to  overall  identified  unique  phosphopeptides  by  each  round  of  Ti02  or 
IMAC  enrichment. 


In  summary,  we  have  determined  that  the  optimal  enrichment  strategy  for 
phosphopeptides  is  to  combine  the  IMAC  and  TiCC  to  enrich  different  subsets  of 
phosphopeptides  from  the  whole  cell  lysate.  We  are  currently  preparing  a  manuscript  on 
these  results  that  will  shortly  be  submitted  to  the  Journal  of  Proteome  Research,  as  a 
follow-up  to  our  manuscript  from  last  year.  The  approach  in  this  part  will  allow  us  to 
include  the  studies  on  phosphoproteomic  changes  caused  by  CDH1  knock  down  as  well, 
which  will  further  enrich  our  conclusions. 
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3.  (Task  Id)  Cellular  responses  with/without  reducing  CDH1  expression.  (8/1/2013- 
9/30/2013,  2/1/2014-5/31/2014,  7/1/2014-current) 

In  order  to  find  the  conditions  where  CDH1  has  clear  effects  on  drug  resistance, 
we  treated  cells  with  three  different  anti-cancer  drugs  including  irinotecan,  oxaliplatin, 
and  5-FU  with  3D  culture  system.  With  irinotecan  treatment,  at  lower  concentration,  the 
CDH1  knock-down  cells  haven’t  shown  much  resistance  to  irinotecan  treatment,  however, 
when  the  concentration  reaches  higher  than  64  pM,  both  Gr  and  B1  cell  lines  that  target 
different  parts  of  CDH1  showed  clearly  higher  cell  viability  after  drug  treatment 
compared  with  control  cells  (Vc)  (Figure  8A).  With  oxaliplatin  treatment,  cells  show 
similar  reaction  as  to  irinotecan  treatment.  At  lower  concentration,  the  CDH1  knock¬ 
down  cells  haven’t  shown  much  resistance  to  oxaliplatin  treatment,  however,  when  the 
concentration  reaches  higher  than  16  pM,  both  Gr  and  B1  cell  lines  showed  clearly  higher 
cell  viability  after  drug  treatment  compared  with  control  cells  (Vc)  (Figure  8B).  In 
comparison  to  irinotecan  and  oxaliplatin,  CDH1  knock  down  cells  do  not  show  any 
resistance  to  5-FU  treatment  at  any  concentration  (Figure  8C). 

The  drug  resistance  test  with  3D  systems  took  longer  than  we  expected  due  to  the 
optimization  of  test  methods.  We  originally  used  agarose-coated  surfaces  for  3D  culture. 
However,  because  the  drugs  can  diffuse  to  the  agarose  layer  that  will  affect  the  actual 
drug  concentration,  the  culture  system  with  agarose-coated  surface  introduced  errors  and 
was  not  stable  between  different  batches  of  tests.  On  the  other  hand,  the  agarose-coated 
surface  cannot  be  directly  used  for  cell  viability  assays  and  the  3D  structures  together 
with  the  cells  that  detached  from  the  3D  structures  need  to  be  transferred  to  new  plates 
for  cell  viability  assays.  This  complication  will  introduce  even  more  errors  to  the  test. 
After  generating  a  substantial  amount  of  data  with  the  agarose-coated  culture  system 
without  being  able  to  identify  any  patterns  of  reaction  to  drug  treatment,  we  concluded 
that  the  agarose-coated  culture  system  is  not  a  very  suitable  system  for  testing  cell 
reactions  to  drug  treatment.  We  then  tried  the  ECM  coated  surface.  On  the  ECM  coated 
surface,  the  cells  also  automatically  grow  into  3D  spheroids,  however,  it  is  hard  to  control 
that  each  well  only  has  a  single  spheroid.  And  also,  even  though  the  ECM  surface  is  a 
thinner  layer,  it  still  has  variability  in  drug  diffusion.  We  finally  switched  to  Ultra-Low 
Attachment  (ULA)  round  bottom  plates,  which  perfectly  solved  all  the  problems  in  our 
tests.  These  plates  do  not  need  extra  coating  with  agarose,  and  can  be  directly  used  for 
cell  viability  assays  without  transferring  cells  into  other  plates.  With  the  ULA  plates,  we 
finally  got  the  above  results  shown  in  Figure  8. 


PAGE  |  18 


Figure  8 
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Fig  8.  Cell  viability  after  (A)  irinotecan,  (B)  oxaliplatin  and  (C)  5-FU  treatment. 
Blue  line  indicates  control  cells  (Vc),  red  (Gr)  and  green  (Bl)  lines  indicate  CDH1 
knock  down  cells. 
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To  further  confirm  the  effects  of  CDH1  knock-down  on  drug  resistance,  we  tried 
other  4  different  siRNAs  that  targets  CDH1.  And  we  tested  2  anti-cancer  drugs  including 
oxaliplatin  and  gefitinib.  As  shown  in  Figure  9,  siRNA-1  and  siRNA-4  clearly  show 
higher  cell  viability  after  drug  treatment,  which  further  confirmed  that  knock-down  of 
CDH1  will  help  cancer  cells  acquire  drug  resistance. 


Figure  9 
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Fig  9.  Cell  viability  after  (A)  oxaliplatin  and  (B)  gifitinib  treatment.  Neg-si 
indicates  cells  transfected  with  siRNA  that  has  non-specific  targets.  Si-1,  Si-2,  Si-3, 
and  Si-4  indicate  cells  transfected  with  different  siRNA  sequences  that  targets 
distinct  parts  of  CDH1  to  induce  CDH1  knock  down. 
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4.  (Task  2b  &  2c)  Collect  mass  spectral  data  for  comparison  of  3D-structures  with 

or  without  CDH1  knockdown  induction  after  drug  treatment  with  the  condition 

confirmed  in  Id,  and  data  analysis  with  MASCOT  search  engine  to  identify  affected 

protein.  (2/1/2014-6/30/2014) 

Currently  we  are  still  in  the  step  of  sample  preparation  and  mass  spectral  data 
collection  stage  for  the  analysis  of  pathways  related  to  CDH1  knockdown  caused  drug 
resistance.  While  optimizing  the  drug  treatment  system  and  conditions,  we  analyzed  the 
differences  between  2D  and  3D  culture  system  as  an  extension  of  the  current  project.  The 
experimental  procedures  are  shown  in  Figure  10.  2D  and  3D  cultured  cells  were 
differently  labeled  with  “light”  or  “heavy”  amino  acids  and  lysed  to  extract  proteins  and 
then  mixed  together.  After  alkylation  and  trypsin  digestion,  the  peptides  were  divided 
into  2  parts.  One  part  was  fractionated  with  high  pH  reverse  phase  column  for  proteomic 
analysis.  The  other  part  was  treated  with  IMAC  to  enrich  phosphopeptides,  and  the 
enriched  phosphopeptides  were  further  fractionated  with  high  pH  reverse  phase  column 
for  phosphoproteomic  analysis  (Figure  10). 


Figure  10 


Fig  10.  Flow  chart  of  experiment  methods  for  proteomic  and  phosphoproteomic 
analysis  to  compare  2D  and  3D  cultured  cells. 
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The  results  are  summarized  in  Figure  11.  We  identified  a  total  of  5834  protein 
groups,  including  2478  phosphoprotein  groups  and  8555  phosphopeptides.  Among  these 
identified  results,  3348  were  unique  protein  groups,  1471  were  unique  phosphoprotein 
groups,  and  6488  were  unique  phosphopeptides.  2468  unique  protein  groups,  979  unique 
phosphoprotein  groups,  and  2069  unique  phosphopeptides  were  identified  in  both 
biological  replicates  (Figure  11A).  Of  all  these  identifications,  703  protein  groups  were 
identified  in  both  non-phosphorylated  and  phosphorylated  states  (Figure  11B). 


Figure  11 
A 

TotalPro 


Distribution  of  Proteome/Phosphoproteome 


UniquePro 

OverlapPro 

QuantPro 

TotalPhos 

UniquePhos 

OverlapPhos 

QuantPhos 


□ 


Fig  11.  Summarization  of  identified  proteomic  and  phosphoproteomic  data.  (A) 
Distribution  of  total  protein  groups,  unique  protein  groups,  phosphoprotein  groups, 
and  unique  phosphoprotein  groups.  (B)  Overlap  between  unique  protein  groups  and 
phosphoprotein  groups. 


The  data  quality  was  further  confirmed  in  Figure  12  and  Figure  13.  Figure  12 
shows  the  data  quality  of  the  mass  spectrometric  proteomic  analysis.  The  correlation  of 
biological  replicate  1  and  biological  replicate  2  showed  a  linear  trend  and  the  R  squared 
value  is  0.816,  indicating  high  correlation  between  the  2  biological  replicates  (Figure 
12A).  The  overlap  between  biological  replicate  1  and  biological  replicate  2  shows  1728 
quantifiable  protein  groups  (Figure  12B),  and  most  of  the  quantifiable  protein  groups  did 
not  show  much  change  in  the  2D  or  3D  cultured  cells  (Figure  12C  &  D).  We  chose  2 
fold  as  a  cutoff  level  to  determine  whether  a  protein  has  different  expression  levels  or  not 
in  the  2D  and  3D  cultured  cells. 

The  data  quality  of  phosphoproteomic  analysis  is  shown  in  Figure  13.  Similar  to 
proteomic  analysis,  the  correlation  of  biological  replicate  1  and  biological  replicate  2 
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showed  a  linear  trend  and  the  R  squared  value  is  0.655  for  phosphopeptides,  indicating 
some  correlation  between  the  2  biological  replicates  (Figure  13A).  The  overlap  between 
biological  replicate  1  and  biological  replicate  2  shows  2069  quantifiable  phosphopeptides 
(Figure  13B),  and  most  of  the  quantifiable  phosphopeptides  didn’t  show  many  changes 
in  the  2D  or  3D  cultured  cells  (Figure  13C  &  D).  We  chose  2  fold  as  the  cutoff  level  to 
determine  whether  a  phosphopeptide  has  different  phosphorylation  level  or  not  in  2D  and 
3D  cultured  cells. 


Figure  12 
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Fig  12.  Data  quality  of  protein  groups.  (A)  Correlation  between  2  biological 
replicates.  (B)  Overlap  between  2  biological  replicates.  (C,  D)  Fold  change 
distributions  of  quantifiable  protein  groups. 
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Figure  13 


Fig  13.  Data  quality  of  phosphoprotein  groups.  (A)  Correlation  between  2 
biological  replicates.  (B)  Overlap  between  2  biological  replicates.  (C,  D)  Fold 
change  distributions  of  quantifiable  phosphoprotein  groups. 


From  the  above  analysis,  we  identified  225  protein  groups  and  136 
phosphopeptides  (representing  78  phosphoproteins)  that  are  up-regulated  in  3D  cultured 
cells,  and  116  protein  groups  and  422  phosphopeptides  (representing  245 
phosphoproteins)  that  are  down-regulated  in  3D  culture  cells.  The  network  analysis 
shows  that  many  of  the  up-regulated  proteins  are  interconnected  (Figure  14A),  and  many 
of  the  down-regulated  proteins  are  interconnected  as  well  (Figure  14B).  More  detailed 
analysis  is  still  in  process. 
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Figure  14 


Fig  14.  Network  analysis  of  up-  or  down-regulated  protein  groups.  (A)  Interactions 
among  proteins  that  are  up-regulated  in  3D  cultured  cells.  (B)  Interactions  among 
proteins  that  are  down-regulated  in  3D  cultured  cells. 


CONCLUSION 

In  the  second  year  of  research,  we  successfully  optimized  the  SILAC  methods  and 
phosphoproteomic  analysis  methods.  The  optimized  method  has  increased  the 
phosphopeptide  identification  for  2  folds.  We  also  optimized  the  system  for  testing  drug 
effects.  Four  different  drugs  including  irinotecan,  oxaliplatin,  5-FU,  and  gefitinib  were 
tested  with  CDH1  knock  down  cells  for  drug  resistance  effects.  Cells  with  down 
regulated  CDH1  levels  show  resistance  to  all  drugs  at  high  concentration  except  5-FU. 
The  mass  spectrum  data  for  drug  resistance  pathways  are  under  collection  and 
investigation.  We  also  analyzed  the  proteomic  and  phosphoproteomic  differences 
between  2D  and  3D  cultured  cells. 

In  conclusion,  during  the  second  year,  we  accomplished  most  of  the  work  we 
proposed  in  our  original  proposal,  and  finished  the  tasks  that  got  delayed  during  the  first 
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year.  The  second  year’s  research  is  a  good  extension  of  the  first  year  research,  and 
provided  us  the  foundation  for  our  future  work. 


IMPACT 

1.  Successfully  optimized  SILAC  conditions. 

2.  Successfully  optimized  methods  for  phosphoproteomic  analysis. 

3.  Successfully  confirmed  the  drug  concentrations  that  CDH1  knock  down  show  clear 
drug  resistance  effects. 

4.  Generated  data  and  successfully  analyzed  the  differences  between  2D  and  3D  cultured 
colorectal  cancer  cells. 

5.  Successfully  started  sample  preparation  and  mass  spectral  data  collection  steps  for  the 
analysis  of  pathways  related  to  CDH1  knockdown  caused  drug  resistance. 


CHANGES/PROBLEMS 

The  main  problem  we  encountered  during  the  second  year  of  the  research  is 
determining  the  drug  treatment  conditions  that  CDH1  knock-down  cells  show  drug 
resistance.  We  originally  used  agarose-coated  surfaces  for  3D  culture.  However,  because 
the  drugs  can  diffuse  to  the  agarose  layer  that  will  affect  the  actual  drug  concentration, 
the  culture  system  with  agarose-coated  surface  introduced  errors  and  was  not  stable 
between  different  batches  of  tests.  On  the  other  hand,  the  agarose-coated  surface  cannot 
be  directly  used  for  cell  viability  assays  and  the  3D  structures  together  with  the  cells  that 
detached  from  the  3D  structures  need  to  be  transferred  to  new  plates  for  cell  viability 
assays.  This  complication  will  introduce  even  more  errors  to  the  test.  After  generating  a 
substantial  amount  of  data  with  the  agarose-coated  culture  system  without  being  able  to 
identify  any  patterns  of  reaction  to  drug  treatment,  we  concluded  that  the  agarose-coated 
culture  system  is  not  a  very  suitable  system  for  testing  cell  reactions  to  drug  treatment. 
We  then  tried  the  ECM  coated  surface.  On  the  ECM  coated  surface,  the  cells  also 
automatically  grow  into  3D  spheroids,  however,  it  is  hard  to  control  that  each  well  only 
has  a  single  spheroid.  And  also,  even  though  the  ECM  surface  is  a  thinner  layer,  it  still 
has  variability  in  drug  diffusion.  We  finally  switched  to  Ultra-Low  Attachment  (ULA) 
round  bottom  plates,  which  perfectly  solved  all  the  problems  in  our  tests.  These  plates  do 
not  need  extra  coating  with  agarose,  and  can  be  directly  used  for  cell  viability  assays 
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without  transferring  cells  into  other  plates.  With  the  ULA  plates,  we  finally  solved  the 
problems  and  successfully  determined  the  drug  treatment  conditions.  However,  this 
caused  some  delay  of  the  second  year  research  and  the  collection  of  proteomic  and 
phosphoproteomic  data  is  still  in  process. 


PRODUCTS 

Manuscripts: 

Yue  XS,  Hummon  AB.  Combining  of  IMAC  and  Ti02  enrichment  methods  to  increase 
phosphoproteomic  identifications,  manuscript  in  preparation. 

Yue  XS,  Hummon  AB.  Proteomic  and  phosphoproteomic  analysis  of  differential  protein 
expression  and  phosphorylation  patterns  in  2D  and  3D  cultured  colorectal  cancer  cells, 
manuscript  in  preparation. 


Presentations: 

Yue  XS,  Hummon  AB.  Proteomic  and  phosphoproteomic  analysis  of  differential  protein 
expression  and  phosphorylation  patterns  in  2D  and  3D  cultured  colorectal  cancer  cells 
(American  Society  for  Mass  Spectrometry  Annual  Meeting;  Baltimore,  MD;  June, 
2014) 

Yue  XS,  Hummon  AB.  Modeling  of  EMT  process  in  3D  cultures  and  proteomic  analysis 
of  differential  protein  expression  in  EMT  related  drug  resistance  in  colorectal  cancer  cells 
(American  Association  for  Cancer  Research  Special  Conference;  June,  2014). 

Yue  XS,  Hummon  AB.  Proteomic  and  phosphoproteomic  analysis  of  differential  protein 
expression  and  phosphorylation  patterns  in  2D  and  3D  cultured  colorectal  cancer  cells 
(Harper  Day  Symposium;  Notre  Dame  IN,  April,  2014) 


PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 

Unchanged 
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SPECIAL  REPORTING  REQUIREMENTS 

Not  applicable. 


APPENDICES  None 
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